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ARTICLE INFO ABSTRACT 


Stimulation of serotonergic neurotransmission by psilocybin has been shown to shift emotional biases away from 22 
negative towards positive stimuli. We have recently shown that reduced amygdala activity during threat pro- 23 
cessing might underlie psilocybin’s effect on emotional processing. However, it is still not known whether psilo- 24 
cybin modulates bottom-up or top-down connectivity within the visual-limbic-prefrontal network underlying 25 
threat processing. We therefore analyzed our previous fMRI data using dynamic causal modeling and used Bayes- 26 
ian model selection to infer how psilocybin modulated effective connectivity within the visual-limbic-prefrontal 27 
network during threat processing. First, both placebo and psilocybin data were best explained by a model in 28 
which threat affect modulated bidirectional connections between the primary visual cortex, amygdala, and later- 29 
al prefrontal cortex. Second, psilocybin decreased the threat-induced modulation of top-down connectivity from 30 
the amygdala to primary visual cortex, speaking to a neural mechanism that might underlie putative shifts to- 31 
wards positive affect states after psilocybin administration. These findings may have important implications 32 
for the treatment of mood and anxiety disorders. 33 

© 2015 Published by Elsevier Inc. This is an open access article under the CC BY-NC-ND license 

( http : //creativecommons.org/licenses/by-nc-nd/4.0/) . 
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1. Introduction 

Serotonin (5-hydroxytryptamine, 5-HT) is an important neuro- 
transmitter within neural networks related to emotion processing. 
We have recently shown that 5-HT2A receptor activation by psilocybin 
(4-phosphoryloxy-N,N-dimethyltryptamine) attenuates amygdala acti- 
vation in response to threat-related visual stimuli in healthy volunteers 
and that the reduction of amygdala blood oxygen level-dependent 
(BOLD) signal is related to psilocybin’s mood-enhancing effect 
(Kraehenmann et al M 2014). Here, we addressed the hypothesis that 
connectivity changes between the amygdala (AMG) and visual and 
prefrontal cortical (PFC) areas contribute to the observed effects of 
psilocybin on threat processing previously observed (Kraehenmann 
et al., 2014). This hypothesis is based on evidence showing that the 
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processing of threat-related visual stimuli may be modulated via 53 
feedback connections from the amygdala to the visual cortex (Furl 54 
et al., 2013). Such top-down input from the amygdala to the visual 55 
cortex may be an important mechanism at the interface between 56 
emotion processing and visual perception — given that the amygdala 57 
has been implicated in tuning visual processing to allocate resources 58 
towards sensory processing of - and coordinating responses to - emo- 59 
tionally salient stimuli (Morris et al., 1998). Furthermore, processing of 60 
threat signals may be modulated via inhibitory feedback connections 61 
from the PFC to the AMG (Hahn et al., 2011; Aznar and Klein, 2013). 62 
Using DCM for fMRI, Sladky et al. (2015) recently analyzed the effects 63 
of the selective serotonin reuptake inhibitor (SSRI) (S)-citalopram on 64 
amygdala-PFC effective connectivity in healthy volunteers. They found 65 
that the PFC exhibited a down-regulatory effect on amygdala activation, 66 
and that this effect was significantly increased by the antidepressant 67 
(S)-citalopram. Importantly, the inhibitory feedback from the PFC to the 68 
AMG has been found to be correlated with 5-HT2A receptor stimulation 69 
(Fisher et al., 2009). Therefore, it is conceivable that the psilocybin- 70 
induced attenuation of amygdala activation (Kraehenmann et al., 2014) 71 


http://dx.doi.org/l 0.101 6/j.nicl.2015.08.009 

2213-1582/© 2015 Published by Elsevier Inc. This is an open access article under the CC BY-NC-ND license (http://creativecommons.Org/licenses/by-nc-nd/4.0/). 
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might be caused by increased inhibitory connectivity from the PFC to the 
AMG. Finally, given the abundance of feed-forward projections from 
visual input regions (e.g. primary visual cortex, VI ) to the AMG (Pessoa 
and Adolphs, 2010) and from the AMG to the PFC (Volman et al., 2013), 
bottom-up connectivity changes may also contribute to psilocybin’s 
effects on threat processing. 

To test these hypotheses, we analyzed the functional magnetic 
resonance imaging (fMRI) data of our previous study (Kraehenmann 
et al., 2014) using dynamic causal modeling (DCM) (Friston et al., 
2003) and Bayesian model selection (BMS) (Stephan et al., 2009). 
DCM is a general framework for inferring hidden mechanisms at the 
neuronal level from measurements of brain activity such as fMRI. Recent 
studies have demonstrated its sensitivity to detect pharmacological ma- 
nipulations in fMRI data (Grefkes et al., 2010; Schmidt et al., 2013b); in 
particular, after serotonergic stimulation (Volman et al., 2013). BMS is 
an essential aspect of DCM studies, as it can be used to test competing 
hypotheses (different DCMs) about the neural mechanisms generating 
data. We applied DCM and BMS to address the following questions: 
First, which is the most likely mechanism underlying threat processing, 
(1) threat-induced modulation of bottom-up connectivity, (2) threat- 
induced modulation of top-down connections, or (3) modulation of 
both bottom-up and top-down connections by threat stimuli. Secondly, 
which of these mechanisms - changes in bottom-up or top-down 
connectivity - contributed to the psilocybin-induced reduction of 
AMG (Kraehenmann et al., 2014) and VI activation (Schmidt et al., 
2013a) in response to threat-related visual stimuli. 

2. Methods 

2.1. Subjects 

In total, 25 healthy, right-handed subjects (16 males, mean age 
24.2 =b 3.42 years) with normal or corrected-to-normal vision were re- 
cruited through advertisements placed in local universities. Subjects 
were screened for DSM-IV mental and personality disorders using the 
Mini-International Neuropsychiatric Interview (Sheehan et al., 1998) 
and the Structured Clinical Interview II (First et al., 1997). Exclusion 
criteria were as follows: pregnancy, left-handedness, poor knowledge 
of the German language, personal or first-degree relatives with history 
of psychiatric disorder, history of alcohol or illicit drug dependence, 
current alcohol abuse or illicit drug use, current use of a medication 
that affects cerebral metabolism or blood flow, cardiovascular disease, 
history of head injury or neurological disorder, magnetic resonance im- 
aging exclusion criteria (including claustrophobia), and previous signif- 
icant adverse reactions to a hallucinogenic drug. Subjects were healthy 
according to medical history, physical examination, routine blood anal- 
ysis, electrocardiography, and urine tests for drug abuse and pregnancy. 
The study was approved by the Cantonal Ethics Committee of Zurich 
(KEK). Written informed consent was obtained from all subjects and 
the study was performed in accordance with the Declaration of Flelsinki. 

2.2. Experimental design 

As previously reported (Kraehenmann et al., 2014), the study design 
was randomized, double-blind, placebo-controlled, cross-over. Subjects 
received either placebo or 0.16 mg/kg oral psilocybin in two separate 
imaging sessions at least 14 days apart. The use of psilocybin was autho- 
rized by the Federal Office of Public Flealth, Federal Department of 
Home Affairs, Bern, Switzerland. Psilocybin and lactose placebo were 
administered in gelatin capsules of identical number and appearance. 
A 0.16-mg/kg dose of psilocybin was selected because it reliably induces 
changes in mood and consciousness, but minimally disrupts behavioral 
task performance and reality testing (Studerus et al., 2011 ). Mood state 
was assessed using the using the Positive and Negative Affect Schedule 
(PANAS) (Watson et al., 1988) and the state portion of the State-Trait 


Anxiety Inventory (STAI) (Spielberger and Gorsuch, 1983) before and 132 
210 min after each drug treatment. The scanning experiment was con- 133 
ducted between 70 and 90 min after drug administration to coincide 134 
with the plateau in the subjective effects of psilocybin (Hasler et al., 135 
2004). Subjects were released about 360 min after drug administration, 136 
after all acute drug effects had completely subsided. 137 

2.3. fMRI paradigm: amygdala reactivity task 138 

Inside the scanner, subjects performed an amygdala reactivity task 139 
comprising alternating blocks of emotional (threat and neutral) picture 140 
discrimination tasks. The picture discrimination task was interspersed 141 
with shape discrimination tasks, which served as baseline tasks and 142 
allowed amygdala responses to return to baseline. 143 

Stimulus material for the amygdala reactivity task was obtained 144 
from the International Affective Picture System (IAPS), a standardized 145 
and broadly validated collection of emotionally evocative pictures 146 
(Lang et al., 2005). Stimulus sets of 48 different pictures were arranged 147 
in picture-triplets on a gray background. The stimulus triplets com- 148 
prised the target picture in the upper center position, and two pictures 149 
as potential matching targets on the left and right sides at the bottom of 150 
the slide. Twenty-four pictures were categorized as threat and 24 as 151 
neutral. The threat pictures were aversive, threat-related pictures such 152 
as attacking animals, aimed weapons, car accidents, and mutilations, 153 
and the neutral pictures depicted activities of daily living, portraits of 154 
humans and animals, and everyday objects. 155 

During the emotional picture discrimination task, subjects were re- 156 
quired to select one of the two IAPS pictures at the bottom of the stim- 157 
ulus triplet that matched the target picture at the top of the triplet. 158 
Selection was indicated by pressing one of two buttons on a magnetic 159 
resonance (MR)-compatible response device with the dominant hand. 160 
A shape discrimination task was performed as a sensorimotor control 161 
and baseline task. This required matching of geometric shapes (circles, 162 
ovals, and rectangles) analogous to the picture discrimination task and 163 
was implemented to control for activation due to non-emotional cogni- 164 
five and visual processing. Both tasks were shown as alternating 24-s 165 
blocks without intermittent pauses. Each block was preceded by a 2-s 166 
instruction (“Match Pictures” or “Match Forms”) and consisted of six 167 
target images that were presented sequentially for a period of 4 s in a 168 
randomized order. The experimental design comprised four repetitions 169 
of the sequence threat shapes neutral shapes, cumulating to a 170 
total duration of 420 s for the complete run. Individual trial durations m 
were not determined by the subjects’ responses, and no feedback was 172 
provided regarding correct or incorrect responses. 173 

2.4. fMRI image acquisition and data analysis 174 

Scanning was performed on a 3 T scanner (Philips Achieva, Best, The 175 
Netherlands) using an echo planar sequence with 2.5 s repetition time, 176 
30 ms echo time, a matrix size of 80 x 80 and 40 slices without inter- 177 
slice gap, providing a resolution of 3 x 3 x 3 mm 3 and a field of view 178 
of 240 x 240 mm 3 . 179 

Data analysis was performed with SPM12b ( http://www.fil.ion.ucl. 180 
ac.uk ). All volumes were realigned to the mean volume, co-registered 181 
to the structural image, normalized to the Montreal Neurological Insti- 182 
tute space using unified segmentation (Ashburner and Friston, 2005) 183 
including re-sampling to 3 x 3 x 3 mm voxels, and spatially smoothed 184 
with an 8-mm full-width at half-maximum Gaussian kernel. First-level 185 
analysis was conducted using a general linear model applied to the 186 
time series, convolved with a canonical hemodynamic response func- 187 
tion (Friston et al., 1994). Serial correlations and low-frequency signal 188 
drift were removed using an autoregressive model and a 128-s 189 
high-pass filter, respectively. Single-subject GLM analysis for the two 190 
sessions (placebo and psilocybin) comprised regressors for threat, neu- 191 
tral pictures, and shapes. These conditions were modeled as box-car 192 
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regressors representing the onset of each block type. Subject-specific 
condition effects for threat minus shapes were computed using 
t-contrasts, producing a contrast image for each subject that was used 
as a summary statistic for second-level (between subject) analyses. 

2.5. Dynamic causal modeling (DCM) 

The current DCM analyses (version 12 with SPM12b) are based on 
the GLM analyses of the fMRI data described above (Kraehenmann 
et al., 2014). In DCM for fMRI, the dynamics of the neural states under- 
lying regional BOLD responses are modeled by a bilinear differential 
equation that describes how the neural states change as a function of 
endogenous interregional connections, modulatory effects on these 
connections, and driving inputs (Friston et al., 2003). The endogenous 
connections represent constant coupling strengths, whereas the 
modulatory effects represent context-specific and additive changes in 
coupling (task-induced alterations in connectivity). The modeled 
neuronal dynamic is then mapped to the measured BOLD signal using 
a hemodynamic forward model (Stephan et al., 2007). We explicitly 
examined how the coupling strengths between VI, AMG, and PFC are 
changed by threat during the AMG reactivity task (modulatory effect). 

2.5.1. Regions of interest and time series extraction 

We selected three regions of interest (ROIs) within a right- 
hemispheric network implicated in visual threat processing, based on: 
(1) previously published conventional SPM analyses of these data 
(Fig. 1 ) (Kraehenmann et al., 2014), (2) previous anatomical and struc- 
tural connectivity studies (Freese and Amaral, 2005), and (3) previous 
DCM studies of threat processing using visual stimuli (Volman et al., 
2013). In DCM for fMRI, a neural network is analyzed in terms of direct- 
ed connectivity changes among selected regions of interest. Regions of 
interest are selected based on both a priori knowledge and hypotheses, 
and on significant task-induced activations. We chose a right- 
hemispheric (subgraph) analysis based on our previous GLM analysis 
of psilocybin effects on threat processing (see Table 1, Fig. 3A and B) 
(Kraehenmann et al., 2014). The rationale for this choice was to ask 
whether the observed psilocybin-induced decrease of right amygdala 
activation in response to threat was mediated by top-down connectivity 
changes from the right prefrontal cortex or by bottom-up connectivity 
changes from the right visual cortex. In addition, we limited our DCM 
analyses to a right-hemispheric network or subgraph in view of statisti- 
cal efficiency: it is common practice to test only a small number of 
regions of interest with DCM. Future DCM studies of psilocybin effects 
on threat processing could include the contralateral homologues of 
the regions investigated here, although our previous GLM analysis did 
not motivate a DCM analysis of the left-hemispheric network. 

The ROIs included: rVl (x = 12, y = — 82, z = —7), rAMG (x = 24, 
y=— l,z=— 13), and the right inferior frontal gyrus within the lateral 
PFC (rLPFC) (x = 54, y = 32, z = 20). The coordinates for the rVl, rAMG 
and rLPFC were based on the contrast of threat pictures minus shapes. 


Regional time series from each subject and session were extracted from 
(10 mm) spherical volumes of interest centered on the suprathreshold 
voxel nearest the group maxima. Time series were summarized with 
the first eigenvariate of voxels above a subject-specific F threshold of 
p < 0.01 (uncorrected) within the anatomical areas, as defined by the 
Pick Atlas toolbox. During time series extraction it may happen that a 
subject does not show activation at the group maximum and that the 
nearest suprathreshold voxel lies outside the anatomical regions. By 
additionally using an anatomical mask, we ensured that time series 
were extracted from within a certain distance of the group maxima 
(10 mm), but were not extracted from a region outside the anatomical 
structure (Dima et al., 2011 ). We could not extract an rLPFC time series 
in two subjects due to lack of individual activations fulfilling both the 
above functional and anatomical criteria. Although it is not necessary to 
preclude subjects who did not show significant activations from the 
DCM analysis, the purpose of DCM is to explain observed activations in 
terms of functional coupling. We therefore restricted our analyses to sub- 
jects who showed significant responses under the assumption that their 
data would provide more efficient estimators of connectivity. 

2.5.2. DCM model space 

First, we specified a three-area base model with bidirectional endog- 
enous connections between VI and AMG and between AMG and LPFC 
(Fig. 2A). VI was selected as the visual input region in our models. All 
visual stimuli were used as inputs. These restrictions allowed us to 
define a small model space. The basic model was then systematically 
varied to provide alternative models of the modulatory effect (induced 
by threat stimuli). The three model variants corresponded to the 
three alternative hypotheses about modulatory effects (bottom-up, 
top-down, or a combination of bottom-up and top-down) and allowed 
us to distinguish between the three hypothesized mechanisms under 
the two treatments (psilocybin, placebo) (Fig. 2B-D). 

2.5.3. Model inference 

Using random-effects BMS in DCM12, we computed expected prob- 
abilities and exceedance probabilities at the group-level to determine 
the most plausible of the three model variants for each drug (psilocybin, 
placebo) separately (Penny et al., 2004). The expected probability of 
each model is the probability that a specific model generated the data 
of a randomly chosen subject, and the exceedance probability of each 
model is the probability that this model is more likely than any other 
of the models tested (Stephan et al., 2009). Bayesian model comparison 
rests solely on the relative evidence for different models (as scored by 
the variational free energy). This evidence comprises the accuracy 
(i.e., percent variance explained) minus the complexity (i.e., degrees 
of freedom used to explain the data). The evidence therefore reflects 
the quality of a model in providing an accurate but parsimonious ac- 
count of the data (and is preferred over conventional accuracy measures 
that may reflect overfitting). Finally, we used random-effects Bayesian 
model averaging (BMA) to compute (subject specific) connectivity 
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Fig. 1. Regional effects from the contrast of threat pictures minus shapes within right lateral prefrontal cortex (rLPFC; z = 20) and right amygdala (rAMG; y = — 1) and from the contrast of 
all pictures (threat of non-threat) minus shapes within the right primary visual cortex (rVl; x =12) across both drug conditions (placebo, psilocybin). SPM{t) overlaid on canonical brain 
slices (thresholded at p < 0.001 uncorrected for visualization). 
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til Table 1 

1 1.2 Dynamic causal modeling parameter estimates. 


1 1 .3 

1 1 .4 

Connection 

Endogenous 


Modulation 


Direct input 


Pla 

Psi 

Pla 

Psi 

Pla 

Psi 

1 1 .5 

VI 

+0.023 ± 0.05 

-0.002 ± 0.01 

— 

— 

+0.011 ± 0.12 

-0.003 ± 0.01 

1 1 .6 

VI -> AMG 

+0.036 ± 0.08 

+0.018 ± 0.05 

+0.027 ± 0.37 

+0.024 ± 0.09 

- 

- 

tl.7 

AMG -+ VI 

-0.028 ± 0.09 

+0.031 ± 0.11 

+ 0.526 ± 1.05 

+ 0.030 ± 0 . 14 * 

- 

- 

1 1 .8 

AMG 

-0.007 ± 0.02 

-0.002 ± 0.01 

— 

— 

- 

- 

1 1 .9 

AMG -> LPFC 

+0.005 ± 0.08 

-0.005 ± 0.06 

+0.103 ± 0.22 

+0.023 ± 0.11 

- 

- 

tl.10 

LPFC -> AMG 

-0.002 ± 0.05 

+0.008 ± 0.00 

-0.394 ± 1.12 

-0.157 ± 0.76 

- 

- 

tl.ll 

LPFC 

-0.014 ± 0.04 

-0.001 ± 0.00 

- 

- 

- 

- 


288 estimates (weighted by their posterior model probability) across all 

289 three models separately for psilocybin and placebo. This conservative 

290 analysis allowed the drug effect to be expressed in all connections and 

291 their threat related modulations, whereby we were able to establish 

292 significant effects in relation to intersubject variability using classical 

293 statistics at the between subject level. 

294 2.5.4. Parameter inference 

295 To evaluate the effect of psilocybin on endogenous connections and 

296 their modulation by threat stimuli, BMA values were entered into two 

297 separate 2-way repeated measures ANOVA with factors drug (psilocy- 

298 bin, placebo) and connection type (endogenous parameters: VI, 

299 VI -+ AMG, AMG VI , AMG, AMG LPFC, LPFC AMG, LPFC; modu- 

300 latory parameters: VI AMG, AMG VI, AMG -» LPFC, LPFC AMG). 

301 Where the ANOVA null hypothesis of equal means was rejected, we 

302 used the post-hoc test (Duncan’s multiple range tests). A paired t test 

303 was further applied to compare direct inputs into VI across both treat- 

304 ments. A p value of less than 0.05 was assumed as statistically 

305 significant. 

306 2.5.5. Correlation with behavioral and mood measures 

307 To investigate correlations between psilocybin-induced changes of 

308 effective connectivity and behavior or mood, the psilocybin-induced 

309 connectivity changes were correlated using Pearson correlations with 

A: base model 



psilocybin-induced changes in behavioral measures (reaction time, 310 
accuracy) and mood scores (PANAS positive affect, PANAS negative 311 
affect, STAI state anxiety). 312 

3. Results 313 

3. 1 . Model inference with Bayesian model selection 314 

Under both psilocybin and placebo, the full model outperformed all 315 
other models with an exceedance probability of 97% (placebo) and 62% 316 
(psilocybin), respectively (Fig. 3). This optimal model comprised bidi- 317 
rectional endogenous connections between VI and AMG, and between 318 
AMG and LPFC, with threat modulating both forward and backward 319 
connections. 320 

3.2. Parameter inference 32 1 

To compare connectivity across drug treatments, the subject-specific 322 
parameter estimates were averaged over the three models for each 323 
treatment using BMA. The endogenous parameters, their threat induced 324 
modulations, and direct inputs from the BMA are shown in Table 1 . Cou- 325 
pling or connectivity in dynamic models is measured in terms of Hz, 326 
where a strong baseline or endogenous connection would typically be 327 
between 0.1 and 0.5 Hz. This means that one can regard the effective 328 
connectivity as a rate-constant. In other words, a strong connection 329 

B: bottom-up model 



C: top-down model 


D: full model 




Fig. 2. Model specification. A, Basic structure of the three-area model: visual stimulus presentation drives VI activity, which is bidirectionally connected to AMG, which in turn is bidirec- 
tionally connected to the LPFC. B, Bottom-up model: the modulatory effect of threat is only mediated via bottom-up connections from VI to AMG to LPFC. C, Top-down model: the mod- 
ulatory effect of threat is only mediated via top-down connections from LPFC to AMG to VI. D, Full model: the modulatory effect of threat is mediated via both bottom-up and top-down 
connections between VI and AMG, and between AMG and LPFC. 
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Placebo 


Psilocybin 


Bayesian Model Selection RFX 



Models 


Bayesian Model SeteOon RFX 



Models 


5 


Bayesian Model Selection RFX 



Models 


Bayesian Model Setedon RFX 



1 2 


Models 


Fig. 3. Results of Bayesian model selection. Bar charts show the expected model probabilities (A, B) and exceedance probabilities (C, D) of the bottom-up model (1 ), the top-down model 
(2), and the full model (3) for the placebo (left) and psilocybin (right) treatment. Notably, the full model with threat-induced modulation of bidirectional connections is the winning model 
for both the placebo and psilocybin treatment. 


causes a large rate of increase in the target region, with respect to 
activity in the source region. The inverse of the connection strength 
can therefore be interpreted in terms of a time constant (i.e., how long 
it would take for a source to increase activity in a target). 

The re was no main effect of drug (Fi >2 2 = 3.10, p = 0.09, r)p = 0.12), 
but a significant main effect of connection type (F 3>66 = 3.94, p = 0.01, 
rfi =0.15), and a significant drug by connection type interaction 

( 1^3,66 = 2.84, p = 0.04, rf v = 0.1 1 ) on modulatory coupling parameters. 
Post-hoc tests on the drug by connection type interaction showed that 
the threat-induced modulation of AMY VI connectivity was signifi- 
cantly reduced after psilocybin compared to placebo administration 
(p = 0.01; Duncan’s multiple range test corrected) (Table 1). There 
was no significant effect of psilocybin on endogenous or input parame- 
ters (all p > 0.05). 

Parameter estimates were obtained from Bayesian Model Averaging 
for placebo (Pla) and psilocybin (Psi), mean =b standard deviation. 
Statistically significant differences between placebo and psilocybin 
treatments (p < 0.05 Duncan corrected for multiple comparison) are 
printed in bold and marked by an asterisk; VI = primary visual cortex; 
AMG = amygdala; LPFC = lateral prefrontal cortex. 


33. Correlation with behavioral and mood measures 350 

We assessed correlations between (psilocybin-placebo) modulatory 351 
coupling changes for the AMG VI connection from BMA and (psilocy- 352 
bin-placebo) changes of behavioral measures (reaction time, accuracy) 353 
and of mood scores (PANAS positive affect, PANAS negative affect, STAI 354 
state anxiety). We found no significant correlations (all p > 0.05). 355 

4. Discussion 356 

In this study, we analyzed the fMRI data of our previous psilocybin 357 
study (Kraehenmann et al., 2014) using DCM, an established framework 358 
enabling tests of directed (effective) connectivity. We were interested 359 
whether psilocybin modulated effective connectivity within a network 360 
implicated in threat processing during an amygdala reactivity task. In 361 
particular, our aim was to differentiate between psilocybin-effects on 362 
bottom-up, top-down, and bidirectional connectivity during threat- 363 
processing within a visual-limbic-prefrontal network. There were two 364 
main findings from our study: Firstly, both placebo and psilocybin 365 
data were best explained by a model in which threat affect modulated 366 
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bidirectional connections between VI, AMG, and LPFC. Secondly, 
psilocybin - compared to placebo - substantially reduced the modulato- 
ry effect of threat on the top-down connection from the AMG to VI . This 
implies that psilocybin attenuates amygdala-dependent top-down 
tuning of visual regions during threat processing. 

Our BMS finding that the full model, which is characterized by bidi- 
rectional modulatory effects of threat on visual-limbic-prefrontal con- 
nectivity, outperformed both the bottom-up and the top-down model, 
is in line with previous DCM studies (Herrington et al., 2011; Goulden 
et al., 2012). In these studies, BMS consistently favored models, which 
implement modulatory effects on both bottom-up and top-down con- 
nections during negative emotion processing. The winning model in 
our study contained reciprocal connections between VI and AMG 
(VI ** AMG) and between AMG and LPFC (AMG ~ LPFC). Both 
VI ^ AMG and AMG ^ LPFC reciprocal connections are critically in- 
volved in negative-emotion processing (Herrington et al., 2011; 
Goulden et al., 2012). In fact, it has been shown that visual threat 
perception may be enhanced through a re-entry mechanism of feed- 
forward connections from VI to AMG and feedback connections from 
the AMG to VI (Herrington et al., 2011). Furthermore, visual threat 
perception may be increased through feed-forward connections from 
the AMG to LPFC (Lu et al., 2012) and attenuated through inhibitory 
feedback connections from the LPFC to AMG (Volman et al., 2013). Al- 
though BMS did not directly compare model fits from different datasets 
(e.g. placebo, psilocybin), our model selection results indicate a consis- 
tent mode of visual threat processing during placebo and psilocybin 
treatments; namely, via modulation of both bottom-up and top-down 
connectivity across the visual-limbic-prefrontal hierarchy. 

Our main finding was that psilocybin (compared to placebo) 
reduced the modulatory effect of visual threat on the top-down connec- 
tion from the AMG to VI. In both humans and animals, visual threat 
poses a strong salience signal, which needs to be processed efficiently 
and therefore binds attentional resources (Pessoa and Adolphs, 2010). 
The “tuning” of visual regions via feedback projections from the AMG 
during threat processing is an important mechanism underlying visual 
threat processing and may enhance perception of visual threat signals 
(Morris et al., 1998). In addition, the AMG has been closely linked to 
salience processing and may, via top-down predictive signals, guide 
bottom-up information processing (Vuilleumier, 2015). Therefore, the 
amygdala may actually determine the affective meaning of visual 
percepts by its effects on sensory pathways — an effect which mainly 
occurs subconsciously and which may be greatly amplified in psycho- 
pathological conditions, such as anxiety disorders or depression. In 
this context, increased AMG reactivity may lead to an increased 
attentional focus on negatively valenced environmental or social stimuli 
and thus effectively blocks out the processing of positive information 
(Disner et al., 2011). This is especially relevant for hallucinogenic 
drugs such as psilocybin, because there has been a close and psycho- 
therapeutically interesting relationship between visual perception and 
affective processes during hallucinogen-induced states (Leuner, 1981). 
The psilocybin-induced attenuation of top-down threat signaling from 
the amygdala to visual cortex may therefore lead to decreased threat 
sensitivity in the visual cortex. This mechanism may crucially underlie 
the previously observed decrease of behavioral and electrophysiological 
responses in the visual cortex to threat stimuli during psilocybin 
administration (Vollenweider and Kometer, 2010; Schmidt et al., 
2013a) and may explain the psilocybin-induced shifts away from 
negative towards positive valence during emotion processing 
(Kometer et al., 2012). In line with the notion that attenuation of the 
top-down connection from the AMG to visual cortex may reduce threat 
processing, a recent study showed that habituation to visual threat 
stimuli may parallel attenuation of top-down connectivity from the 
AMG to visual cortex (Herrington et al., 2011). In addition, it has been 
found that hyper-connectivity between the AMG and visual cortex 
may underlie increased threat processing and anxiety (Frick et al., 
2013). 


Given the relevance of LPFC in regulating AMG activity during threat 433 
processing, and given previous studies showing that serotonergic stim- 434 
ulation may increase inhibitory top-down connectivity from LPFC to 435 
AMG (Pessoa and Adolphs, 2010; Volman etal., 2013), we hypothesized 436 
that psilocybin-induced reduction in AMY activity might be due to an 437 
increased LPFC AMG top-down connectivity during threat processing. 438 
However, psilocybin did not appear to increase top-down connectivity 439 
from LPFC to AMG in the current analysis. Two reasons might account 440 
for this. First, the source of the psilocybin-induced reduction of AMG 441 
activity, as observed in our previous GLM analysis (Kraehenmann 442 
et al., 2014), might not reflect an increased top-down effect from 443 
LPFC, but rather a suppression of recurrent interactions with visual 444 
areas mediated by a reduced top-down connectivity with the visual 445 
cortex. The synaptic basis of this reduced top-down modulation might 446 
reflect a direct effect of psilocybin in the amygdala: amygdala neurons 447 
abundantly express 5-HT2A receptors, and DOI and other 5-HT2A ago- 448 
nists produce direct effects in the amygdala (Rainnie, 1999). In addition, 449 
a directly decreased AMG reactivity would result in a reduced load on 450 
the LPFC to regulate AMG activation. This view is supported by a recent 451 
DCM study showing that increased AMG-related load on the PFC yields 452 
subsequent responses in the PFC to regulate the AMG (Volman et al., 453 
2013). Second, the AMG might be regulated by prefrontal cortical re- 454 
gions other than the LPFC, such as the medial PFC (MPFC), the anterior 455 
cingulate cortex (ACC), or the orbitofrontal cortex (OFC), which have 456 
also been related to the ‘aversive amplification’ circuit (Robinson et al., 457 
2013). For example, Sladky et al. (2015) recently analyzed the effects 458 
of the selective serotonin reuptake inhibitor (SSRI) (S)-citalopram on 459 
amygdala-OFC effective connectivity in healthy volunteers. They 460 
found that the OFC exhibited a down-regulatory effect on amygdala 461 
activation, and that this effect was significantly increased by the antide- 462 
pressant (S)-citalopram. Although Sladky et al. used a similar threat- 463 
inducing amygdala reactivity task (Hariri et al., 2002) and likewise 464 
tested the effects in healthy volunteers, their study procedures differ 465 
substantially from our study, both in terms of task design (e.g. face stim- 466 
uli instead of pictures, scrambled control stimuli, longer baseline condi- 467 
tions) and in terms of drug administration (e.g. chronic and repeated 468 
instead of acute and single treatment). Therefore, it is not easy to disam- 469 
biguate task- from drug-specific effects in terms of PFC involvement and 470 
our DCM might have missed top-down effects from PFC on the AMG. 471 
However, given the cognitive task requirements in our task - where 472 
subjects were not explicitly required to evaluate or regulate their emo- 473 
tional responses to the threat stimuli - and given that the GLM analyses 474 
(Kraehenmann et al., 2014) did not show significant BOLD responses in 475 
the MPFC, ACC, or OFC, one might argue that top-down effects from 476 
other prefrontal regions are unlikely. Overall, both the hallucinogen 477 
psilocybin and the non-hallucinogen (S)-citalopram may normalize 478 
amygdala hyper-reactivity to threat-related stimuli; leading to their 479 
antidepressant and anxiolytic efficacy, but psilocybin appears to 480 
regulate emotion processing and mood by acting on network 481 
interactions which are different from classical antidepressants such as 482 
(S)-citalopram, such as the affective regulation of visual information 483 
processing shown here. 484 

4. 1 . Limitations and future directions 485 

There are some limitations to be considered in the present study. We 486 
used a fairly simplistic neuronal network underlying threat related 487 
effective connectivity. There are also other brain regions involved in 488 
threat processing, such as the ACC, the OFC, or the fusiform gyrus 489 
(Robinson et al., 2013), but that we did not include in our present 490 
model for reasons of parsimony and based on our a priori hypotheses. 491 
Furthermore, to maximize statistical efficiency, we only considered 492 
right-hemispheric networks in our DCM analyses. Therefore, top- 493 
down connectivity from the left LPFC to the right AMG might have 494 
been missed. Given the importance of the left LPFC in regulating the 495 
right AMG during emotion processing and in serotonergic modulation 496 
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(Outhred et al., 2013), we cannot exclude this possibility. Therefore, fur- 
ther effective connectivity studies using tasks that differentially recruit 
left and right prefrontal cortical regions during threat processing, are 
needed. 

4.2. Conclusion 

This effective connectivity study shows that a decrease of top-down 
connectivity from the AMG to the visual cortex underlies the psilocybin 
effect on visual threat processing. This result suggests that decreased 
threat sensitivity in the visual cortex during emotion processing may 
explain the potential of psilocybin to acutely shift emotional biases 
away from negative towards positive valence: the capacity of the visual 
cortex to process multiple stimuli is limited and hence top-down sup- 
pression of negative stimuli enhances the processing of positive stimuli 
(Kastner et al., 1998). This may have important therapeutic implications 
for mood and anxiety disorders, where over-loading with negative 
stimuli and persistence of negative cognitive biases is a central feature 
(Disner et al., 2011). In post-traumatic stress disorder, for example, 
psilocybin might help inhibit fear-responses during exposure-based 
psychotherapy, which might facilitate therapeutic progress. 
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